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We present a dynamical lattice calculation with 2 flavours for bottomonium states with an 
additional gluonic excitation. Using improved actions for the quarks and gauge flelds at a lattice 
y—^ • spacing of a « 0.1 fm, we find 10.977(61) (62) GeV for the energy of the lowest lying bbg-hjhrid, 

' where the first error is statistical and the second denotes the systematic uncertainty due to the 

. determination of scale. In a parallel quenched simulation we demonstrate explicitly that vacuum 

04 ' polarisation effects are less than 10% of the splitting with the ground state. 

PACS: 11.15.Ha, 12.38.Gc, 12.39.Hg, 12.39.Jh, 12.39.Mk, 14.40.Nd 

' The intense experimental search for particles with exotic quantum numbers has triggered theoretical interest in 
CO . hadrons which contain an excitation of the gluon flux [Q. The lattice approach to QCD has already resulted in 
consistent predictions for the masses of hybrid quarkonia from first principles (2|-||]. The recent results suggest 
that the lowest lying hybrids are 400-500 MeV above the phenomenologically interesting BB threshold. However, for 
numerical simplicity, all those calculations were performed in the quenched approximation, where vacuum polarisation 
effects have been ignored. In this paper we present our results from a full calculation in which dynamical sea quarks 
' are included. Dynamical results for light hybrids have been reported by the SESAM Collaboration ||^ . Here our main 
CO , goal is to establish whether the dynamical treatment of the light sea quarks and the heavy valence quarks will result 
in large shifts of the hybrid levels. 

On the lattice, hybrid mesons are particularly difficult to treat as their excitation energies are much larger than for 
conventional states. This makes it very difficult to resolve the exponential decay of the meson propagator on coarse 
, lattices. It has been demonstrated that such problems can be circumvented on lattices with a finer resolution in the 
7^ ' temporal direction but the implementation of dynamical sea quarks on anisotropic lattices has yet to be achieved. 
Q-i. This work is part of a larger project to study sea quark effects in QCD on isotropic lattices For the gluon sector 
^ ' we chose a renormalization-group improved action, which is written in terms of standard plaquettes, TrP^^, and 1x2 
rectangles, TvR^,^: 
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^ = 4 {coTrP^. + ciTri?^,} (cq + 8ci = 1) . (1) 



Our prescription ci = —0.331 is motivated by an RG-analysis of the pure gauge theory For the light sea quarks 
in the gauge field background we chose a lattice formulation which removes both the doublers and 0{a) discretisation 
errors: 

Sq^q {/^ + rrig) q + a q I^^ q - Cs^{a) a ^q (J^^^F^^ q . (2) 

Further details and results for the conventional light hadron spectrum can be found in ||^. 

Here we study heavy hybrid states as they are of particular relevance for ongoing experiments at B-meson factories. 
To this end we implemented a non-relativistic approach (NRQCD) for the heavy b quarks on a fine lattice with 
a « 0.10 fm. This approach is well suited, owing to the small velocity of the quarks within the flat hybrid potential. 
Unlike calculations for the spin structure in quarkonia we expect only small corrections for spin-avcragcd quantities 
from higher order relativistic terms. This assumption has already been tested explicitly for charmonium Q|. To 
leading order we describe the forward evolution of the heavy quark as a diffusive problem 
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The parameter n was introduced to stabilise the evolution against high-momentum modes. We chose n = 2 throughout 
this analysis. As it is common practise, we have also included the correction terms C7 and Cg to render the evolution 
equation (||) accurate up to C(a^), classically. Radiative corrections induce terms of 0{aa'^) and we applied a mean- 
field improvement technique to reduce such errors as first suggested by Lepage et al. We decided to divide all 
gauge links by the average link in Landau gauge: uq — (0|(l/3) [/p|0)L, which is a suitable choice to reduce the 
unphysical tadpole contributions in lattice perturbation theory. 

Non-relativistic meson operators can be constructed using the standard gauge-invariant definitions of [Q. For the 
magnetic hybrid this amounts to an insertion of the colour-magnetic field, Bi, into the bilinear of two-spinors: 

3 

'H{x)^Yl '/'^(^) ■ (4) 

1=1 

In the following we will denote as IB and 2B the ground state and the first excitation onto which this lattice operator 
can project. We also employ fuzzed link variables and several different smearings for the quark fields. This results in 
different projections of the source operators onto the ground state and it is useful when extracting higher excitation 
with the same quantum number. 

For our study we calculated the bottomonium spectrum on 400 dynamical configurations at (/?, k) — (2.10, 0.1382) 
on a 24^ x 48 lattice. This corresponds to the lightest sea quark mass of our full data set at (3 — 2.1 and we measured 
niTr/mp = 0.5735(48) For the comparative quenched analysis we accumulated 192 independent configurations at 
(3 = 2.528. This coupling was chosen so as to match the lattice spacing of the dynamical run. In both cases we find 
a ~ 0.11 fm from the string tension y/a — 440 MeV. Here we take the IP — IS splitting to set the scale. As expected, 
such a definition results in slightly different lattice spacings for our two data sets. 

Finite size effects are known to be small for heavy quarkonia and they have been explicitly checked for IS, IP and 
IB in charmonium |^. Since here we study the bottomonium system on even larger lattices {L « 2.5 fm), we do not 
expect any volume dependence for the ground state masses in our analysis. The simulation parameters for our two 
data sets are collected in Table |. 

The results for hadron masses are obtained from correlated multi-exponential fits to different smearings, a, and 
timeslices, t: 

"fit 

C(a,i)-(Mt(t)Af„(0))-^a,(a)e-"'* . (5) 

1=2 

In Table || we present the results from multi-exponential fits (ngt < 4). 

As it is standard in lattice calculations with heavy quarkonia, we have tuned the bare quark mass in Equation 

so as to reproduce the experimental value of M^in/ilP — IS) = 21.5, where aMkin is determined for the non- 
relativistic dispersion relation of the S-state on the lattice. In Figure [l] we demonstrate the mass independence of 
the spin-independent ground states at (/3 = 2.528, Nf — 0), along with their higher excitations. This independence is 
important as it allows us to extract predictions at a slightly non-physical point. It also means that our result is stable 
against possible radiative corrections to the quark mass. Since in NRQCD calculations all energies are measured with 
respect to the ground state, IS", we introduced the ratio Rx — {X ~ IS) / {IP — IS) to quote the normalised splitting 
with respect to the IP — IS. 

Here we focus on the ratio Rib = {IB — IS) /{IP — IS) which determines the mass of the lowest lying hybrid. 
The velocity expansion of NRQCD works very well for the hybrid states, owing to the slow quarks in a fiat hybrid 
potential. Therefore we do not expect significant changes for mass of the spin-averaged hybrid state due to higher order 
relativistic corrections. To leading order in the NRQCD Hamiltonian, the singlet state of Equation^ is degenerate with 
the corresponding spin-triplet states (0 ^, 1 ^,2 ^). These arc all the states with zero orbital angular momentum 
and they include the exotic combination 1 due to the coupling of the spin to the giuon angular momentum. A 
near degeneracy was also reported for hybrid states with additional orbital angular momentum |^ . Such states give 
rise to more exotic states, such as 2^ . The above-mentioned degeneracies will be lifted once higher order relativistic 
corrections are re-introduced into the NRQCD Hamiltonian [T^JT^ . 

The main uncertainty in lattice calculations of hybrid excitations so far is the absence of dynamical sea quarks. 
Previous estimates of quenching errors frequently referred to the uncertainty in the determination of the lattice spacing 
as a limiting factor of quenched simulations. However, it is not clear a priori, whether the high gluon content of the 
hybrid itself would cause large shifts in its mass once dynamical sea quarks are introduced. Comparing the results in 
Table |l[ we find, perhaps surprisingly, that this is not the case. Indeed, our quenched estimate for Rib = 3.43(45) is 
in excellent agreement with the dynamical simulation. Rib — 3.59(14). This means that quenching errors are smaller 
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than our statistical errors of about 10% for this quantity. This is a pleasing feature for lattice calculations and it is 
also confirmed by another dynamical calculation which could not resolve any change of the static hybrid potential. 
In this picture one can also understand why relativistic corrections and discretisation errors are so small - in the flat 
potential the quarks move very slowly and are widely separated, hence they are not very sensitive to the details of 
the lattice cutoff. This is in contrast to the S-state splittings, where dynamical effects are to be expected because 
of their sensitivity to the physics at short distance scales. Here we observe a 5% shift of R2S, in line with an earlier 
calculation at similar lattice spacing | p5[ . It is important to realize that this accounts for only part of the dynamical 
effects in the real world. As it has been argued before our sea quark mass is sufficiently small for the purpose 
of bottomonium calculations, but it is apparent that large gluon momenta can also excite strange quarks from the 
vacuum - a full description ought to include three dynamical flavours. 

For a short-range quantity like i?2S it is equally important to study scaling violations. Our previous work shows 
that discretisation errors may also account for a fraction of the remaining < 10% deviation from experiment. We leave 
a more precise determination of those two different sources of systematic errors to future studies. For Rib we have 
convincing evidence from previous simulations |^,^ that at a « 0.1 fm scaling violations are negligible for spatially 
large hybrid states. In the following we take the deviation of R2S from its experimental value as a conservative 
estimate for residual systematic errors in our 2-flavour simulation. 

In Figure || we plot our results against Nf, along with previous quenched estimates from different lattice spacings 
and a variety of isotropic/ anisotropic lattice actions. The stability of all these results and the good agreement with 
our value for 2-flavour full QCD provides support to the arguments above. Converting our dynamical result into 
dimensionful units we quote 10.977(61)(62) GeV for the lowest lying hybrid, where the first error is statistical and 
the second denotes the systematic uncertainty in the determination of the lattice spacing. 

Taking the conservative average of all points in Figure ||, we find 11.02(18) GeV as the lattice prediction for the 
mass of the lowest spin-averaged bottomonium hybrid. We want to stress the importance of this result, as it shows 
that even after introduction of two dynamical flavours the lowest lying hybrid will be above the BB threshold (10.56 
GeV), where a number of experiments are currently running. In fact our prediction is intriguingly close to the BB** 
threshold (« 11.0 GeV), below which hybrid states are thought to be very narrow. 

Our results show that the error from the quenched approximation is small within the statistical error of isotropic 
lattice calculations (w 1% of the bound state energy). To study systematic errors for hybrid states at an even 
higher level of accuracy is possible on anisotropic lattices, but it remains to be seen how dynamical quarks can be 
implemented in such a framework. 

The calculations were done using the supcrcomputing facilities at the Center for Computational Physics at the 
University of Tsukuba. This work is supported in part by the Grants-in-Aid of Ministry of Education (No. 09304029). 
TM, HPS and AAK are supported by the JSPS Research for Future Program, and SE and KN are JSPS Research 
Fellows. 
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(Ns,Nt) 


ao- [fm] 


aip-is [fm] 


{aniQ, n) 




config. 


(2.10, 0.1382) 


(24,48) 


0.11293(44) 


0.10403(80) 


/o 0/t 0\ 


0.854604 


400 


(2.528,quenched) 


(24,48) 


0.11319(97) 


0.0883(16) 


/I no o^ /o 0/1 o^ 


0.876700 


192 


TABLE I. Simulation parameters. For the dynamical run 


we measured am-^/amp 


= 0.5735(48). 


i0,Nf) 




(2.528, 0) 


(2.528, 0) 


(2.10, 2) 




experiment 


amq 




1.92 


2.24 


2.24 










4.008(16) 


4.628(17) 


4.642(18) 






IP- IS 




0.2057(34) 


0.2007(37) 






439.37(13) MeV 


Mkin [GcV] 




8.74(15) 


10.35(19) 


8.805(76) 




9.46030(26) 


(2S-1S)/(1P-1S) 




1.461(51) 


1.459(48) 


1.389(22) 




1.2802(16) 


(3S-1S)/(1P-1S) 




3.12(27) 


3.15(23) 


3.255(98) 




2.0353(24) 


(2P-1S)/(1P-1S) 




2.28(18) 


2.21(19) 


2.336(60) 




1.8186(20) 


(3P-1S)/(1P-1S) 




4.36(40) 


4.47(42) 


4.94(17) 






(1B-1S)/(1P-1S) 




3.43(45) 


3.63(39) 


3.59(14) 






(2B-1S)/(1P-1S) 




5.8(1.1) 


6.1(1.1) 


6.83(52) 







TABLE II. Normalised splittings with respect to IP— IS. The experimental values are taken from rC^Si and the spin-averaged 
nP. The dynamical results come from our lightest sea quark mass (k = 0.1382). 



4 



8.0 

7.0 

6.0 

w 5.0 

^ 4.0 
w 

>< 3.0 
2.0 
1.0 
0.0 



• X=2S 
CX=3S 
■ X=2P 
□ X=3P 

♦ X=1B 
<>X=2B 



<> 



i 



1 



i 



1 



8.5 



9.5 
M.n[GeV] 



10 



10.5 



11 



FIG. 1. Mass independence of spin-avcragod bottomonium spectrum. Here we plot the excitations against the kinetic mass 
of ^Si from two different quark masses. The measured energies axe 8.74(15) GeV and 10.35(19) GeV. The hybrid results are 
shifted for clarity and the experimental Upsilon mass is shown as a vertical dashed line. 
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FIG. 2. Dependence of Rx on numbers of dynamical flavours, Nf. The experimental values are plotted at Nf = 3. For Rib 
we also show quenched results from other groups as triangles (offset for clarity). 
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